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The conductivity in quasi two-dimensional systems is calculated using the quantum kinetic equa- 
tion. Linearizing the Lenard-Balescu collision integral with the extension to include external field 
dependences allows one to calculate the conductivity with diagrams beyond the GW approximation 
including maximally crossed lines. Consequently the weak localization correction as an interference 
effect appears here from the field dependence of the collision integral (the latter dependence some- 
times called intra-collisional field effect). It is shown that this weak localization correction has the 
same origin as the Debye-Onsager relaxation effect in plasma physics. The approximation is applied 
to a system of quasi two-dimensional electrons in hetero-junctions which interact with charged and 
neutral impurities and the low temperature correction to the conductivity is calculated analytically, 
ft turns out that the dynamical screening due to charged impurities leads to a linear temperature de- 
pendence, while the scattering from neutral impurities leads to the usual Fermi-liquid behavior. By 
considering an appropriate mass action law to determine the ratio of charged to neutral impurities 
we can describe the experimental metal-insulator transition at low temperatures as a Mott-Hubbard 
transition. 



I. INTRODUCTION 

The low temperature conductivity of quasi two- 
dimensional systems like MOSFET structures or hetero- 
junctions reveals a surprising metal to insulator 
transition 1 . This critical review of the theoretical ap- 
proaches has concluded that this phenomenon is insuffi- 
ciently explained. The experimental data shows a pro- 
nounced transition from insulating behavior at low densi- 
ties to a metallic behavior at high densities. The generic 
feature of the metal-insulator transition is the rapid 
change from insulating to conducting behavior when the 
density is increased very slightly at low temperatures. 
This density driven metal-insulator transitions are usu- 
ally referred to as Mott transitions 2 . The characteris- 
tic feature of the Mott-Hubbbard transition is that an 
increase in the effective mass is directly responsible for 
increasing resistivity p — m/e 2 nT while the Anderson 
scenario would assume a vanishing relaxation time nr. 
Measurements of the effective mass 3 seem to support a 
Mott-Hubbard transition rather than the Anderson tran- 
sition quantitatively explained in Ref. 4 . 

In this paper we want to return to the original idea of 
the Mott transition in that a bound state is resolved with 
increasing density due to pressure ionization. We will 
show that a quantitative description of the experimental 
results can be achieved if one calculates the interplay of 
weak localization and trapping due to charged impurities 
as well as the scattering with neutral impurities on the 
same theoretical footing. 

Weak localization as a quantum interference effect 
is intensively studied in the literature 5-9 . The maxi- 
mally crossed diagrams lead to a diffusive pole which al- 
lows one to extract weak localization corrections to the 
conductivity 8,10-12 (and citations therein). This paper 



is devoted to an alternative route to weak localization. 
We will linearize the quantum kinetic equation derived 
at a lower level of approximation but with external fields 
in order to create higher order diagrams. The main aim 
is to show that weak localization has the same origin as 
the interference effect known from the Debye -Onsager 
relaxation effect in plasma physics. 

In the following part of the introduction we will out- 
line the model we want to use. This will give a summary 
about the kinetic approach adopted in this paper. The 
many-body approximation to be applied in this paper 
is also specified and it is clarified how higher order dia- 
grams are generated by linearization, which is presented 
in detail in appendix A. In the second section we calcu- 
late the electrical conductivity from the kinetic equation 
approach and present results for the relaxation time and 
the relaxation function. This relaxation function sum- 
marizes the quantum corrections to the conductivity due 
to interference effects. The Lenard-Balescu kinetic equa- 
tion used as a starting point in this section is derived 
in appendix B. The third section discusses the result- 
ing conductivity formulae and shows that the relaxation 
function is identical to the one calculated for the weak- 
localization corrections. For comparison with experiment 
we include scattering from neutral impurities besides the 
scattering with charged impurities which are worked out 
in detail in appendix C. The fourth section summarizes 
the results and gives an outlook. In appendix D we dis- 
cuss the polarization function in quasi two-dimensional 
systems. The additional appendices E and F present cal- 
culations of the integrals used during the paper. 
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A. Outline of the model 

As a model, we assume a quasi two-dimensional 
Coulomb potential, where the field lines are three dimen- 
sional but the motions of the particles are restricted to 
two dimensions. The effective potential takes then the 

form V a b = 4:ne a e b h/ ^Jq% + g 2 where for instance e a = e 

for electrons and e b — —Ze for charged ions. Our ap- 
proach is conveniently based on the kinetic equation for 
the one-particle distribution f a (k,t) normalized to the 
area density 



dk 



(27rft) ; 



:/a(M) 



n„ 



(1) 



where the spin degeneracy is denoted by s. The kinetic 
equation for that distribution function is 



d t fa + e a EV k /o 



b 



(2) 



where correlations are covered by the corresponding col- 
lision integral X. This collision integral is explicitly field 
dependent due to the distortion of two-particle correla- 
tions. We consider the conductivity in a system described 
by a local equilibrium distribution 



/«(M) = e 



(k- Pa (f)) 2 - 



+ 1 



(3) 



with the mean mass-motion of the charged particle p a (t). 
The center of mass motion is at rest, which means that 
the total sum of currents J2 Jb = E ""m^" = with 

b b 

density n& and mass mb- In the following we will restrict 
ourselves to a two-component system. The generalization 
towards multicomponent systems is straightforward. 

From the collision integral (2) we have two sources 
of linear response: A term proportional to the current 
j = — e a n a p a /m a and a term proportional to the field. 
By multiplying (2) with k and integrating, the balance 
equation for the momentum reads 

/dk 

tiSE -i 
= -n Q e a E— + n a p a r 

hi 

such that the current balance takes the form 



<9tja - 



i-fl.-r-'j. 



(4) 



(5) 



(6) 



According to the distortion of the Fermi function (3) we 
have a linear response 



and the stationary conductivity j = crE is 



SE 
E 



fa - fa = -Pa - <9 k J fo 



(7) 



and we can represent the current relaxation time and 
the relaxation effect for a scattering with impurities (C4) 
with the potential V s as 



SEj 
E 



'wS(e k - tk-q) 

x 1 ft v 

2 efe— ek-q 



(8) 



The relaxation time r is coming from the term propor- 
tional to the current in the collision integral. This is due 
to on-shcll scattering represented by the delta function in 
(8). Besides this the conductivity becomes renormalized 
by the explicit field dependence of the collision integral, 
which is an interference effect and corresponds to 
an off-shell scattering as expressed by the derivative of 
the principal value in (8). The latter effect has been 
the subject of various investigations for nondegenerate 
plasmas and is known as Debye-Onsager relaxation ef- 
fect. The field dependence and quantum form is dis- 
cussed in 13 . For a plasma system this Debye- Onsager 
relaxation effect 14-20 was first derived within the theory 
of electrolytes 21-25 . * 

The theoretical calculations of conductivity in reduced 
dimensions is the topic of intensive investigations. These 
concern rigid two dimensional electron systems 27 and 
quasi- two dimensional systems 28 . In the latter study a 
three-dimensional system was considered, where the par- 
ticles can only scatter in two dimensions leading to a 
cylindrical Fermi surface. The Born approximation and 
contact interaction result in a resistivity which has a lead- 
ing low temperature behavior as — cx aT 2 (l + 61nT). 
We will show that the Coulomb interaction with the dy- 
namical screening results in a linear order as the leading 
term. This has been repeatedly reported in the liter- 
ature both from an experimental and theoretical point 
of view. Numerical calculations of Coulomb scattering 
rates from impurities predict a linear temperature depen- 
dence of the mobility in silicon inversion layers 29 ' 30 . This 



* Debye derived a limiting law of electrical conductivity 1 
which stated that the external electric field E on a single 
charge Z = 1 is diminished in an electrolyte solution by the 
amount SE/E — where e is the elementary charge, E the 
electric field strength, T is the temperature of the plasma and 
k is the inverse screening radius of the screening cloud. This 
law is interpreted as a deceleration force which is caused by 
the deformed screening cloud surrounding the charge. Later 
it has been shown by Onsager 22 that this result has to be cor- 
rected to SE/E — 3(2 ^y2 )T if the dynamics of ions (Z = 1) 
is considered. The linear response theory reproduce this On- 
sager result 17 ' 18,20 . The kinetic theory leads to the Onsager 
result if asymmetric screening 13 is applied while the symmet- 
ric treatment leads to the Debye result 19 ' 26 ' 20 . 
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was attributed to the collisional level broadening in the 
screening function. Related results have been obtained in 
Ref. 31 where a significant suppression of the temperature 
dependence of the screening function was found. An ana- 
lytical investigation of screening in quasi two-dimensional 
systems was given in Ref. 32 where a linear temperature 
term in the conductivity was reported. 

In this paper we want to investigate the effect of 
Coulomb screening on the conductivity. We will derive 
exact analytical results which show that due to dynami- 
cal screening the leading low temperature contribution to 
the conductivity is linear. In contrast, the static screen- 
ing leads to a quadratic temperature dependence typical 
for the Fermi liquid. 



turn, 5 can be expressed as a variation of the selfenergy 
£ with respect to the Green's function (bare line) G. For 
the dynamical screened approximation used here we give 
the corresponding results in figure 1. One sees that in 
principle maximally crossed diagrams are accounted for. 

Instead of the variation of the nonequilibrium Green's 
function we can use a proper reduction of the latter ones 
towards a kinetic equation. The linear response obtained 
from this kinetic equation including all external field ef- 
fects is accounting then for higher order diagrams in a 
convenient way. As such, we obtain weak localization 
effects by proper linearization of the collision integral. 

II. KINETIC EQUATION AND CONDUCTIVITY 



B. Many body approximation used in this paper 

For the calculation of the conductivity we want to ac- 
count for quantum interference effects such like weak lo- 
calization. This means we have to include maximally 
crossed diagrams at least 5 ' 7 ' 12 . Besides the direct cal- 
culation of diagrams, we can employ the philosophy of 
variational techniques. In Ref. 33 it was described how 
one can use the variation of nonequilibrium Green's func- 
tions with respect to an auxiliary external field to create 
higher order diagrams in the response function. This 
makes use of the known variation technique summarized 
in appendix A. 



X = i 




Ji + 



= — 

8 G 
= + 



= + 




FIG. 1. The density response function (first line) in terms 
of the particle- hole T- matrix. The latter one can be expressed 
in the second line as a sum of irreducible graphs which are 
given in terms of the selfenergy. The latter one is used in 
screened potential V = V + VJ1V (last line) approximation. 

The response function can be given in terms of the par- 
ticle hole T-matrix which in turn can be represented by 
the sum of irreducible graphs, 5, according to figure 1. In 



The kinetic equation corresponding to the dynamical 
screened approximation is the quantum Lenard - Balescu 
equation, which has been derived for high fields in 13,34 . 
A sketch of the derivation is in appendix B 



dt Ja 



eE ^/« = C(M-/° ut (iM). (9) 



The collision-in integral is [I out is given by / <-» 1 — / 
and L< <-> L>] 



OO 



du> 
2^ 



X cos 



t e a Eqr 2 



h 2m a h 

x/ (k-q-e Er,t-T)(l-/ a (k-e Er,t-r)) 
n<(q,c,t-lr) 
\£(q,u,t- ±t)| 

with the free density fluctuation (B4) 



(10) 



n<(q,w,i) = -2 



dp 



dr 



X cos 



b t e b Eqr 2 



(hw-e b p + e b p+q )- + 



h 2nibh 



1 1 
x/b(p + q,i- 2 T )(! _ fb(P,t- 2 T ))- 



(11) 



Here we use the sum over species explicitly. 

The nested form of the Lenard - Balescu collision inte- 
gral (10) is computationally advantageous. It tells that 
the collision of particle a with momentum k on a parti- 
cle b with momentum p into (a, k — q) and (b,p + q) can 
be represented equivalently as a collision of the particle 
a with k on a hole a with k — q by a dynamic plasmon 
emission which is considered as a particle-hole fluctua- 
tion of particles b. For static screening £(q, 0, t), equation 
(10) reduces to the kinetic equation for statically screened 
Coulomb potentials in high electric fields (CI) 35 ' 34 . 
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Now we calculate the frequency integral in (9) and 
(10) analytically using the identity (El). We summarize 
the result of the frequency integration in the momentum 
dependent function W(q) which takes the explicit form 
(E10) for quasi two-dimensional systems. 

Performing the balance equation for the current (4) we 
obtain in linear response the interference term or relax- 
ation function 



6E 
~E 



OO 

x J dTT 2 I s {a,T)(-hi<{q,Q)) - (a<-> 6) 



(12) 



and the relaxation time as 



with 



n, 



T -^=^(R(a,b)- e ^R(b,a)) (13) 

OO 

x J dTTl s (a,T)(-^U<(q,0)). 



(14) 



Here the static free density fluctuation reads 



n>(g,o) = n<( g ,o) 



' ,q f p (l-f p )S(e p -e p+q ) 



(15) 



and the integral 

I s (a,r) = 



rlk 



x sin (e fc _ 9 - e fe ) 



(16) 



will be calculated in appendix F. 

We remark that the interference effect (12) vanishes 
for identical scattering partners, e.g. electron-electron 
correlations. 

When calculating the explicit form of the relaxation 
time and relaxation function we employ charge neutrality 
e a n a + ebnb = and restrict ourselves to the case of single 
charge ions e b = —e. The case of higher charged ions is 
also available but is more involved. The generalization 
to systems with additional particles species is straight 
forward. 



1. Relaxation time by charged impurities 

Introducing the dimensionless integration variable q = 
ZpfaU we obtain the relaxation time to lowest order in 
temperature 



1 = E 



m a e a pf a s b 
irh 4 n a 
l 

dy- 



K a Ti 



e b n b 



yn a 



£ + 1 J (V + Ka) ijl^y 1 



(17) 



where we have introduced the abbreviation f = m \ e \ . 

m a e a 

Since the momentum integration y ^ 1 is restricted by 
the low temperature expansions to values below 2pf a 
the low temperature expansion of the inverse screening 
length (Ell) becomes a constant 



h 2 



l 



e&Sfc 



(18) 



We distinguish here between the temperature of electrons 
T a and the temperature of the ions T b which could mime 
nonequilibrium effects. 

The integral in (17) can be easily calculated 



R = t- 1 



m a e„ 



8s b f T a , T b 




k' 



K a K a 



2 yi— 7^ (i + yr^)^ 

where we will use the abbreviation 



2p fa 2Hpfa 



roe/ 



(19) 



(20) 



witllTn 1 



s 2 /8/h 3 during the paper. Further- 



more, we distinguish in (20) between the inverse screen- 
ing length where ions are included, n' a , and where they 
are neglected n a 



< = «a(l + V^~) 



(21) 



We expand the above result for large and small K a 
which corresponds also to the small and large density 
limits 



"-fa tfaj l+£ 

1 _ * + 2 + r K >-3) 
^n a + K 2 a \nf + o(K 3 a ) 



(22) 
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It is interesting to investigate the limit of large ion 
masses, £ — > oo and T b = 0, which would correspond to 
the charged impurity limit. One gets from (22) for large 

R=^^n a . (23) 

If we compare this with the neutral electron-impurity 
scattering result (CIO), we see significant differences. 
While the statically screened result shows a Fermi-liquid 
behavior of const and T 2 terms, the dynamically screened 
result (23) leads to a linear temperature dependence. 



The low density (large K a ) expansion as well as the high 
density (small K a ) expansion read ' 

SE _ _2sb (T^,Tb_\ _J_ 

j + (In (2<)- l)^+ (O 
1 -5 + 2<-^k? + («?) 

(29) 



III. DISCUSSION 



2. Relaxation function by charged impurities 

Interference effects from the relaxation function (12) 
can be calculated analogously 

~E = 2*p fa hn a (ebTa CaTb) J+l 



y{y + <) 



TZrA Vy^&(y- i±v) (24) 



The small n has been introduced to perform the princi- 
pal value integration V according to (F9). It should be 
noted that the functions of the denominator and nu- 
merator cancel exactly and no restriction on y-integration 
remains. 

We have now to carefully consider the structure 

oo 

f d 2 

z± = J d yfyg^2 i9 v ®{y - 1 ± vn = 



x 1™ J dyfyg'y + (f y g' y - f v g v )\ y= i Tn (25) 



with f y = K,'Jy/(y + n a ) and g y = yjy 2 - 1. Performing 
the integral one sees that the divergent contribution at 
n — > is cancelled exactly by the fg' — f'g term. We 
obtain 



For further progress we use the expansions for large 
K a or low densities. Collecting equations (22) and (29) 
we obtain the conductivity in an analogous form to the 
Bloch-Gruneiscn formulae 



e 2 l-¥ 



(n) 



R{n) 



(30) 



where the dynamic parts comes from the scattering from 
charged impurities (22), (29) 

f = f^ + ^) (l-(ln(2 Ka )-l) \ + o { n-J)) . 

(31) 

We find that both the relaxation time as well as the re- 
laxation function have a linear temperature dependence 
for temperatures lower than the Fermi energy This is in 
agreement with the experimental and theoretical works 
mentioned in the introduction. 



A. Relation to weak localization 

The low density or weak potential limit (30) is now 
interesting to discuss. Using (20) we can rewrite (30) 
into a net relaxation effect [s a = s b = 2, T a = T b ] 



T± = 



1 - k'- 



± 



7r arcoshft' 



«2 _ 1)3/2- 



(26) 



The principal value in (17) is calculated from (Z + +Z_)/2 
and we obtain (with charge e b — — e a ) 



8E 



2s b (T a T b \ £ 
■ — H ) K, 



e fa e/o 



arcosh/i' 



\-n' 2 « 2 -l)^ 



(27) 



^It is clear that for the impurity limit with infinite masses 
£ — > oo we have from (21) and (29) 



5E_ 
~E 







(28) 



in agreement with the physical picture that if the ions cannot 
move the screening cloud cannot deform during the motion 
of the electrons and cannot induce a relaxation effect. This 
is different if the charged impurities do not contribute to the 
screening, K a = K ! a in (21), and we obtain a finite result. 
This case is anticipated here since the neutral static relaxation 
function already lead to finite results (C14). 
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a_ _ _ SE 

t h ( h \ 

where we have introduced the temperature dependent 
relaxation time 1/r = efTro/H 2 = 8efhT/me 4 s 2 . 
This is precisely the localization correction to the 
conductivity 7,6 . Therefore we understand now the phys- 
ical meaning of the relaxation function at low tempera- 
tures. Please note that it vanishes here for small temper- 
atures in contrast to the static result (C14). 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 
T[K] 



FIG. 2. The conductivity versus temperature according to 
(36) as dashed lines. The experimental curves are the solid 
lines 36 . From top to bottom they correspond to densities of 
n a = 6.85, 7.17, 7.25, 7.57, 7.85 x 10 10 cm- 2 . 



B. Comparison with experiments 

In order to describe realistic experiments, we extend 
our model with scattering from static neutral impurities 
(treated in appendix C14). In addition to (30) we have 
static contributions coming from the scattering from neu- 
tral impurities (C10),(C14) as 



R, 



2 3/2 t 



2r 



7T TT 2 T 2 TT 3 

2m 24 e|~16^ 



To(«- 4 ) 



6E n a \2r 



4 - 3 In 2nj _ 4 
1 + -a + o(K i ) 



(33) 



We can safely use the large n a limit since typical den- 
sities of 36 are 7 x 10 10 cm~ 2 and we have 



284.9 



y/n a /7 x 10 



-w cm -2 m 



(34) 



which is a large parameter. 

Now we have two unknown fit parameters in the the- 
ory. These are the energy level Eb — E c — Ed of the 
impurities determining the ratio of impurity density to 
electron densities 37 ' 4 



2irh 2 
m a T 



(35) 



and the ratio of the scattering length to the scattering 
range of neutral impurities C = {r^j From (30) we 

can write the final conductivity formula in the form 



- = 1 - ±e E ^ T 
<J 2T 



= -r 



e 2 ( C 



h VV2T 



T 


T 


2e/a + 


2 £faK 2 a 


T 













111- 



(36) 



where we have used the spin degeneracy of the heavy 
impurities s& = 1 and the temperature X& = 0. The best 
fits to the experimental results 36 are plotted in figure 2. 



We see in figure 2 a clear insulator to metal transition 
for low temperatures when the density is increased very 
slightly. The fitting formula (36) works quite well at all 
experimental densities for low temperatures, however the 
formula fails for higher temperatures. This is because we 
used the low temperature Sommerfeld expansion and the 
Fermi energy is 1.9 K xn a /7 x 10 10 cto~ 2 in this case, 
such that at 0.8K we expect deviations from the leading 
low temperature behavior. Despite this imperfect agree- 
ment with the data, it is quite satisfying that the metal 
to insulating transition can be described completely by 
the scattering with charged and neutral donor impurities 
supplemented by a mass action law. This strongly favors 
the Mott-Hubbard transition picture. 



Fit [K] 




1.02 1.04 1.06 1.08 1.1 

n [7x 10 10 cm- 2 ] 

FIG. 3. The fit parameter Et and £ from (36) versus den- 
sity. 

The best fit parameter are shown in figure 3. We 
see that the effective binding energy of the electrons to 
the impurities as well as the scattering strength decrease 
rapidly when passing through the critical density. This 
is in agreement with the picture that we have pressure 
ionization i.e. a crossing between continuum and bound 
state levels. Therefore, this underlines the Hubbard- 
Mott transition which we have seen already in the in- 
crease of the effective mass 3 as described in 4 . 

A remark concerning the formula (36) should be made 
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here. Instead of the trapping at charged impurities re- 
sulting in the mass action law (35), one could assume in 
principle any other trapping mechanism. In particular, 
we have shown that three-particle bound states can de- 
scribe the experimental data as well 38 . The underlying 
conductivity formula is precisely (36) where the num- 
ber of three-particle bound states diminishes the charged 
impurity density correspondingly. So the basic transport 
mechanisms outlined in this paper here remain the same. 
The only mechanism that is not possible to extract so far 
is that of actual trapping. 
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APPENDIX A: VARIATIONAL TECHNIQUE OF 
LINEAR RESPONSE 



IV. SUMMARY AND OUTLOOK 

Linearizing the Lenard - Balcscu collision integral in- 
cluding all external field dependences allows one to de- 
rive a conductivity on the level of an infinite series of 
diagrams including maximally crossed lines. The field 
dependence of the collision integral yields an interference 
effect which is shown to describe just the weak local- 
ization corrections. This has the same formal origin as 
the Debye-Onsager relaxation effect in plasma and elec- 
trolyte systems. 

For the low temperature regime, it was possible to cal- 
culate the conductivity analytically. It was found that 
the dynamical screening by charged impurities leads to a 
linear temperature dependence of the conductivity, while 
neutral impurities give rise to the usual Fermi liquid be- 
havior. This finding is general for any scattering of light 
particles from heavy particles. Therefore it might also be 
of use for scattering rates in high T c cuprates 39 . 

The comparison with experiment is performed assum- 
ing an appropriate mass action law between the charged 
donor impurities and the neutral ones which are consid- 
ered to be captured electrons. The experimental metal- 
insulator transition can be described quantitatively by 
fitting the effective binding energy and the unknown scat- 
tering strength. 

As noted, a similar quality of description of the ex- 
perimental data was achieved assuming a three particle 
clustering instead of trapping of electrons. The latter 
process, however, relies on the same transport picture as 
outlined in this paper. Only the composition of neutral 
and charged impurities is determined differently. So far 
we cannot determine which process is actually happen- 
ing. In order to achieve this we must study the magnetic 
field dependence, which shows quite unique and remark- 
able features in the experiment. This is left for further 
work. 

To conclude, we suggest that the metal insulator tran- 
sition found in experiments can be described within a 
Mott-Hubbard transition scenario in agreement with the 
effective mass measurements 4 . 



Assuming besides the interaction potential Vw also a 
coupling of an external potential Uw where the numbers 
sign cumulative indices like space time,... coordinates, 
we can express the two particle Green function Gm'2' = 
l/i 2 < Ta^a^a^ > by a variation of the one particle 
Green function G12 = l/i < Ta\a^ > with respect to 
the external potential 40,41 as 



Gl21'2' = G11/G22' T 



6G, 



5U 2 



(Al) 



where the upper sign denotes the Fermi and the lower 
the Bose functions. Using the Dyson equation 



G- 1 = Ga 1 - E - U 



(A2) 



we can calculate the derivative in (Al) and with the help 
of the chain rule and 5G = —G5G~ 1 G, one can express 
the fluctuation function as 



'121'2' 



121'2' 



G\yG<xii 



G 

=FGi2'G 2 i' T G13 crr 34 G41' 



— TGl2'G21' + G 



'5U 2 >i 



34 



13 



5G 



^5262 



G41'. 



56 



(A3) 



Double occurring indices are understood as integrated 
over. With the definition of the occurring vertex function 
we can express this graphically 




- + 




_ 8 XL _ ^ 

~~ 34 ^ m t-t 



5G 



M 

^3645 



56 

FIG. 4. 
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Sometimes it is of advantage to express this density fluc- 
tuation function by the T-matrix. Defining 




T 



FIG. 5. 



we can express 



^12 = + i 



= +i 




+ i 



1 



FIG. 8. 



where we have used the definition of the T-matrix in the 
u-channcl figure 5. The figures (4), (5) and (8) consti- 
tutes the approximation used in this paper (1). 



- + 




FIG. 6. 



APPENDIX B: DYNAMICAL SCREENED 
APPROXIMATION IN NONEQUILIBRIUM 

The selfenergy is given in terms of the dynamical po- 
tential V according to figure 1 



With the help of (Al), the density response to an ex- 
ternal potential can be expressed in terms of the density 
fluctuation function L of (A3). Therefore we remark that 
the density is given by iGn+ =< a^ai >=< h\ >= n\ 
and we have from (Al) and (A3) for the response func- 
tion x 



£<(M,t') = 



dq 



:V<(q,M')G<(k-q,M') 



(Bl) 



where the dynamical potential is expressed within 
Coulomb potentials V a b(q) 



X12 



Sni 

su~ 



22 



= ±i < («i - ni)(n 2 - n 2 ) > . (A4) y< (q, £, t') = ]T V ad (q)C< (q, t, t')V ca (q) 



(B2) 



The last identity follows from the definition of L and un- 
derlines the names density fluctuation function. We see 
now that the linear density variation due to an external 
potential can be expressed as 



-n\ = TL 121+2 +U 2 = -Xi2^22 



(A5) 



where the upper index indicates the order of external field 
dependence. Graphically we can express it as 



X 2 



in 1 

i M 1 




- + 



FIG. 7. 

where we will design the external field as a dotted line 
ending with a cross. According to (A5) and figure 7 we 
can express the first order response function as 



dc 



via the density-density fluctuation 



C<(q,t,t') = 5 at 



dtdt 



x (rr 1 (q,a)n< ( q ,t,l) (ry 1 ( q J,t'). 



(B3) 



Here II is the free density fluctuation or polarization func- 
tion 



n<(q,M') 



dp 



2 G<(p,t,t')G>(p-q,t',t) 



(B4) 



and E r l a the retarded/ advanced dielectric function 

r>( q , t, t') = s(t - 1') ± iG[±(t - ]T v bb (<i) 

b 

x(n>(q,t,t')-n < (q,t,t'))- (B5) 

One easily convince oneself that this set of equations (BI- 
BS) is gauge invariant. The correlation or Green's func- 
tion can be related to the Wigner distribution f a by the 
GKB ansatz 42 

G<(k,r,*)=exp{4(e fc T + ^/)} 



x/ ( k-^!il,t-ii! 

J 1 2 2 



(B6) 



and analogously for G > by replacing / <-» (1 — /). 
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With the help of the gauge invariant formulation of 
Green's function, we can write the kinetic equation for 
the Wigner function /(p, t) = G<(p, R, t,r = 0) finally 35 



d_ 

di 



t-to 



/(k,t)+eEV k /(k,t)= J dr 



r,t--),E<(k-— T,-T,t--) 



G<(k-^r,r,t-I),E>(k-^r,-r,*-I) 



(B7) 



This kinetic equation is exact in time convolutions. This 
is necessary because gradient expansions in time are con- 
nected with linearization in electric fields and conse- 
quently fail 43 . The gradient approximation in space has 
been applied assuming slow varying processes in space 
and we have dropped all R-dependence for simplicity. In- 
troducing (Bl) into the equation for the Wigner function 
(B7) one obtains the kinetic equation (9) with the explicit 
form of collision integral (10). 



APPENDIX C: STATICALLY SCREENED OR 
FINITE RANGE IMPURITY SCATTERING 

Using the static approximation for the dielectric func- 
tion £(q, 0, t) in (10), the kinetic equation for stati- 
cally screened Coulomb potentials in high electric fields 
appears 35 ' 44 

d d 
;fa + eV—f a = Y,Iab 



df 



Iab ^ f^^ks iMb -K-K) 



ft 2 J (2irh) 4 

OG 

V 2 (k a -k' a ,t) J dr cos |(e + e b - e' a - e' 6 )^ 



Er 2 ( e a k a e b k b e a k' a e b k h 



2ft \ m a mb m a m,b 

X {f a >Ml - /a)(l - fb) - fa.f b (l ~ fa>){l ~ fb')} 

(CI) 

describing the scattering of particles a (electrons) with 
other species b with the distribution function fb = f b (k b — 
e b Er 7 T — r). The potential turns out to be the static 
Debye one 



V s (q) 



2ne a e b h 



(q + Kk) 

with the static screening length n given by 
K = ^ 2ire 2 c d^n c 



(C2) 



(C3) 



and the chemical potential /i. 

We will now use this statically screened result in order 
to describe the scattering neutral impurities if we use the 
range of potential r = 1/k and replace the charges by 
the scattering strength g a b = e a e b . 

The calculation of the impurity scattering in quasi-two 
dimensions is now analogously to the Brooks-Hearing re- 
sult for three dimension and starts from the Born col- 
lision integral (CI) which takes for infinite heavy ions 
[m b /m a ->]oo] 

i ab (k a ) = |f j ^^ s(K-K)v s 2 (K-K,t) 

oo 

/, j ( ,,t Er 2 fe a k a e a k' a \\ 

x dr cos < (e a — e„)- — - > 

J \ K a, h 2h \m a m a J J 

o 

Xfb{fa' ~ fa}- 

(C4) 

We assume parabolic bands e = k 2 /2m. 

The relaxation function SE/E would correspond to lin- 
earization of the cos function in (C4) with respect to the 
field while the relaxation time is obtained taking into ac- 
count the linearization with respect to the momentum p a 
of the displaced distributions (3). Cross terms like p a E 
are already of second order response. The result can be 
written in the concise form (8). 



a. Relaxation time by neutral impurity scattering 

In the following we give an explicit calculation. Em- 
ploying the Yukawa or Debye potential (C2) one obtains 
for the relaxation part 

oo 

T Rn \ df 2irm a 2 f cos (a - <p) - cos (f> 

labika) = -Pa-R — 9 ab \ do.- 



(2/c a | sinf | + Hk) 2 



(C5) 



with the angle between k a and the field direction p a de- 
noted by (j). 

The current relaxation time is now obtained by 



naV^ei 1 = j ( 2 ^ 2 k a-fafe(fca) (C6) 

from which one gets [k p = h/2rop] 

In 



-i _ ™ a g 2 ab n b s b f 

-¥J^—J dpdpfo 



n 2 -l ( K 2 -l) 3 / 2 



(C7) 



This current relaxation time in the low temperature 
Sommerfeld-expansion leads to [«j = ft/2r p/ a ] 
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-1 = m a9ab n bSb 

Tt 2 3 / 2 ft 3 n„ 



7T 2 T 2 

+ — — d 2 

+ 24 e 2 K - 



In 



+ 



K 2 - 1) 3/2 
ln /WEM 



where we have used (F9) for the last line. Employing the 
same regularization due to the principal value, (25), we 
end up with 



K -1 



(«2 - 1)3/2 



) 

(C8) 



5E__n_ 
~E ~ n, 



H ( mg ab r Q \ 2 J 2k|+k^ 3«:g In ^k 2 - 1) { 
Q V ft 2 J \(« 2 -l) 2+ («2-l) 8 / 2 /• 

(C13) 



The T 2 dependent term could be given analogously. The 
needed low density or short range expansion reads now 



The second temperature correction is negative and di- 
minishes the positive first part. As long as T < ef the 
net relaxation time is positive and continously falling to 
zero for /c, — > oo. 

The high density expansion reads 



WmwX 2 | 1+ 4-31n2. o | +0(K _ 4) ^ 



E ~ n, 



Ri = r,r 



m a 


e 2 


n a e 2 a 


/l 


*{- 


1 - 


TT 2 


T 2 


24 


e 2 



2 3/; 



s 6 n fc 



(^) 2 .f(l + ^) 2 



2 



APPENDIX D: POLARIZATION FUNCTION IN 

2D 

Here we discuss the properties of low temperature po- 
larization function 



|_Z T ( 2 + (19 + 121n ^ ? ) +o(K 3 ) n K «) = /^ 



d P M 2m 2m j 

ft) 2 e-^-^o 



(Dl) 



(C9) 

and the low density expansion or short range expansion 
2 3 / 2 s b n b /m a .g ah ro^ 2 



The imaginary part is easily rewritten as 



Ri 



{ 1 



h 2 

TT 7T 2 T 2 TT 3 



k 2 2k\ 24 e 2 16k 3 



«f(i + Ve) 2 
+ («r 4 ) 



/dpp jp 1 
(2W /0( ^ } 



2tt 



(CIO) 



v / 9 2 \ X ,P-Q , <f \ 

x < o( w) — o( h U)j 

\ m 2m ; v m 2m 1 

DO 

(q/2 + mLoh/q) 2 



m^ 2 Jde ( 



fe. Relaxation effect by neutral impurity scattering 

The relaxation function is now obtained if we linearize 
(C4) with respect to the external field. We obtain 



(g/2 - muh/g) 2 
- fo{e + 2^ } 



2m 



(D2) 



The energy shifts in the distribution function we absorb 
into an effective chemical potential which should be pos- 
itive in order to obtain nonzero contribution at Sommcr- 
feld expansion 



H<1,E) J, 



X cos drr sin 



q 2 + 2kq 

2mh 



(Cll) 



Meff = £/ - " q -^~ > 



2m 



(D3) 



from which one gets 

oo 

SE n, 



The low temperature Sommerfeld expansion reads than 
f (e) = n((e-^ B )/T) 



dqq 3 V 2 (q) J drr 2 I s 
- /' dy 



E 4irh 5 mn a 

o 

(C12) 



rfe e - /i c ff 

:7l 



T ) = 2 y dxn(x)d x ^Tx + fi cS 

Hetl/T 

oo 

= 2 y — n{x))y/Tx + /i c jf 



Meff/T 



10 



= 2VM^F i - 



2 1 



7T 2 T 2 

24/i cfr 

7t 2 t 2 a 

12 de f 



'fJ-eS- 



(D4) 



Therefore it is enough to know the zero temperature re- 
sult since the T 2 correction are given simply by deriva- 
tives with respect to e/. 






2 4 6 8 10 

FIG. 9. The real (solid) and imaginary (dashed) part of the 
polarization 47r/t 2 /m x II versus x for x = 0.5 (above), x = 1 
(middle) and x = 2 (below) according to (D5). 

Introducing dimensionless coordinates as in 28 
q huo 



2p f 



x 



4c/ 



(D5) 



we get finally for the imaginary part of the polarization 
function 



ImIT(cj, q) 



m 



4irh 2 x 



x je(z - |x + .t 2 |)^/i- (* + ^ 

-e(a: - |z - ^1)^1 -(x-^y 



(D6) 



which is of course the result given in 28 . The correspond- 
ing real part is given by the Hilbert transform according 
to (Dl) 



ReII(w, q) = -2 J 
Using the integral 



duo' Imn(w', q) 
2n lo — lo' 



(D7) 



/(«) 



a — z 



1 > \a\ 



a — sgn(a)Va 2 — T 



(D8) 



we obtain 
Ren = -- 



m 



{/(f +*)-/(^ -*)}■ (D9) 



47T 2 ?l 2 X 

The real and imaginary part is plotted in figure 9. 



APPENDIX E: INTEGRALS OVER DIELECTRIC 
FUNCTIONS 

In order to perform the frequency integration in (9) we 
use a very useful relation, which has been given in 45 



/ 



2tt lo yHl ' 2 



1 



£(q,0)J ' 



(El) 

For the integration of (9) we set H(u>) = 
cos (wt + A) cos (lot 1 + B)lo /lm£ with A and B are the 
remaining content of the cosine functions of (10) and (11). 

Lets first prove the relation (El). We consider the 
following integral including the dielectric function 

1= I ^^Mlmf-V) 



/ 



2ir lo 
dto ( 1 



+ - 



47ri \lo + i-q lo — it] 



H{w){f--n 



(E2) 



where /+ = 1 — 1/ £ and /~ = (/ + )*. In the following we 
will assume that the function H(to) is analytical. Since 
f ± (to) has no poles in the lower/upper half plane and 
vanishes with <~ lo~ 2 for large lo we have the identity 

'^LH{ U )J^L = Tf ± {0)H{0) (E3) 



/ 



2wi 



(lo ± irj) 

and all other combinations of signs vanish. With the help 
of the relation (E3) we compute easily for (E2) 



(E4) 



which proves relation (El). 
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1. Regularization of integration 

We are now going to give explicit forms including the 
dielectric function £ = 1 — V(q)H(q, 10) where the polar- 
ization function II was discussed in the foregoing chapter. 
The forms appearing throughout the paper are 



du) h{uo) 
2^\£(u)\ 2 
duj w/i(w)/Im£(w) 
2ir ui 



Im£-V) 



(E5) 



such that we can apply (E4) with H(u>) = 
wh(u))/lm£(u)). In the case where Im£ = appears 
an ambiguity which we have to remove. We add to the 
particle-hole fluctuation which forms the polarization an 
infinitesimal small classical process 



SlmU = r]u>e 



rj(uj + o(lu 3 )). 



(E6) 



This will make the imaginary part of the polarization 
nonzero everywhere and will not introduce additional 
poles in the upper half plane. Therefore we can apply 
the integration (E4). The corresponding real part ac- 
cording to (E4) will be oc r\ and drop out in the final 
form (E4). 



2. Specific forms 

Now we write down the required forms for (E4). 
One gets the static result from (D6) and (D9) at low 



temperatures 



28 



q < 2p f 



with the chemical potential fi and from (D2) 

Imll(cj, q) = 
m 2 u)hO(2p f - q) ( 8ir 2 T 2 m 2 ' 



(E7) 



irh 2 q^jAp 2 - q 2 
+o(lu 2 ,T 4 ) 



(E8) 



The region where Imll ^ correspond exactly to the 
upper case of (E7). Using this expansion we obtain 



H(0) = lim 



Im£ K 



hq 2 



E 



(E9) 



where we have to keep in mind that the above procedure 
in calculating the frequency integral works only for finite 



Imll or finite particle-hole fluctuations. According to 
(D2) this restricts the later g-integration to values smaller 
than 2pfb respectively. 
We get finally for (El) 



■ 1 ,. hui 



1 



2ite 2 a e 2 tfK{q) 1 



q + hn(q) 



(E10) 



where the screening length is from (E7) 



1 



q < 2p fb 



q > tyfb 



(Ell) 



and since n = j^zp} + o(e C/ / T ) one has d^n = f^p-. 

APPENDIX F: LOW TEMPERATURE 
EXPANSION OF INTEGRALS 

The integrals occurring in (12), (14) and (16) will now 

oo 

be calculated. Using J dr cos xt = ttS(x) we can write 
o 

for (B4) and (15) 

~U<(q,w = 0) = J -^f b {p){l-f b {p + q)) 

oo 

x J dr' cos (e p+q - e p ) ^- 



dp 



2 fb(p){l - fb(p))nh8 (e p+q - e p ) 



{2ttK) 



"' 1 dpMp)(i-f b (p)) 1 ,lr 



Airhq 



m 2 T b j ^ n{x)(l-n(x)) 
nhq J X J%m b (T b x + e fb ) - q 2 



(Fl) 



where we used x — p 2 /2m b T b — €f b /T b and n(x) = l/(c x + 
1) and e/b is the Fermi energy. The last integral is only 
non-zero for negative lower integration limits implying 
tyfb > Q- Expanding the square root in terms of the 
temperature T b we obtain finally 

~n<(q,u = 0) = ~ll>(q,u = 0) = 



12 



mlT b e(2 Pfb g) ( | 8^ 2 m 2 T 2 



(F2) 



The integrals (12), (14) can be tremendously simplified 
observing that (16) can be written 

I s (a, t) = J dkf k (l - fk- q ) sin (e fe _ 9 - e fc ) ^ 
= J dk(fk- q - fk)g ek _ q - ek sin (e fe _ 9 - e k ) ^ 

= J dkfk (g, k+q ~, k sin (e k - e k+q ) I 

-5e fe - £fc _,sin (e fe _ g - e fe ) I) 
= J dkf k sin (e fc - e k - q ) T - {ge k _ q -e k + 9e k -e k _ g ) 
= J dkf k sin(e k - q - e k ) ^ (F3) 



where ^ = l/{e x / T — 1) and we have used k — > — fc 
transformation in coming from the third to the fourth 
equality. 

The time integral (14) can be now represented as a 
derivative of a S function with respect to k. A partial 
integration leads than to 



oo 

/ 



drrl s {a 1 r) 



oo 

/ (2^ /fe / 



drr sin (e k - q - e k ) 



h 



m a f dk 1 1 £ /, \ \ f dx 

W J T d ^ k)) 

q/2 



X\/\ — X 2 



OO 

= f dkd k (kf a (k)) 

q/2 



fc 2 - £ 



3/2 



4V2tt(7 



OO 

/ 

-o/T a 



da:n(x)(l — n(x)) 



e(2 P/a - g) 
Vr a z + a 



(F4) 



with a = e/ — q 2 /8m a . Expanding the argument in 
terms of T a one arrives analogously to (F2) at the result 



oo 

J ' dTTl s (a,r) 



1 



8^mjT 2 

n 2 _ „2A2 



X0(2p /o -g). (F5) 
The integral (12) requires some more care. We use the 

oo 

principal value identity / dr sin tx — V/x and write 

o 



drT 2 I s (a, t) = 



oc 

/ 



^ OO 27T 



COS( 



(F6) 



Using 



2tt 



2tt 



cos</> Va 2 - 1 



0(a- 1) 



(F7) 



we can write after one partial integration 



oc 

/ 



drT 2 I s (a, t) 



ixq ? 



q 2 /8m a T a -e fa /T a 
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/ 



dxn(x)(l — n(x)) 



-Cfa/T 



xJ? 2 - 8ra T a x - 4pj a . 



(F8) 



Expanding again the argument in terms of T a we arrive 
at 



oo 

J nq* 



(F9) 
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